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I.  Introduction 


Recently  active  experiments  from  rockets  and  the  Space  Shuttle  have  been  conducted  to 
study  beam  propagation,  spacecraft  charging,  instabilities,  wave  ('mission,  and  other  space 
plasma  phenomena.  There  experiments  indicate  that  artificial  electron  beams  with  energy  of 
several  ke\  and  currents  about  100  m A  produce  strong  plasma  heating  and  electromagnetic 
radiation  in  the  whistler  frequency  range.  Although  beam-plasma  interactions  for  a  finite- 
radius  beam  system  tire  well  known  for  electrostatic  waves,  electromagnetic  instabilities  of 
a  finite-radius  beam  system  have  not  been  extensively  studied  because  of  their  complexity. 

In  this  study  we  focused  on  the  electromagnetic  instabilities  of  a  finite-radius  electron 
beam  in  a  neutralizing  background.  This  report  gives  the  derivation  of  the  electromagnetic 
dispersion  equation  for  the  finite-radius  beam  system  and  documents  the  numerical  code 
developed  for  solving  the  dispersion  equation.  Appendix  A  presents  a  study  of  electrostatic 
instabilities  of  a  finite-radius  beam  system. 
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II.  Background 


Electron  beam  injection  experiments  have  been  conducted  from  rockets  and  the  Space 
Shuttle  to  study  beam  propagation,  instabilities,  wave  emission,  and  other  space  plasma 
phenomena  in  the  ionosphere  [Winkler,  I960;  Grandel,  1982;  Shawhan  et.  al.,  1984;  Obayashi 
et  al. ,  1984],  Experiments  have  shown  that  electron  beams  can  generate  intense  waves  of 
a  broad  frequency  range.  For  example,  during  the  Spacelab  1  mission,  electron  beams, 
which  were  fired  from  the  Shuttle  with  a  beam  energy  up  to  8  keY  and  a  beam  current  up 
to  300  mA,  produced  broadband  emissions  at  frequencies  less  than  the  electron  cyclotron 
frequency  (js  <  fie)  with  both  electric  and  magnetic  components  [Brghin  et  al.,  1984;  Taylor 
et  al.,  1985).  During  the  Spacelab  2  mission,  the  Plasma  Diagnostics  Package  (PDP),  flying 
around  the  Shuttle  at  distances  of  up  to  300  meters,  again  detected  broadband  emissions 
from  fie  at  about  1  Mz  down  to  30  kHz  while  the  shuttle  emitted  an  electron  beam  with  a 
beam  energy  of  1  keV  and  current  of  50  mA  [Gurnett  et  al.,  1986].  With  both  electric  and 
magnetic  components,  the  broadband  emissions  had  a  clear  funnel  shape  in  the  frequency¬ 
time  spectrogram  similar  to  the  auroral  hiss  reported  by  Gurnett  et  al.  [1983],  Using 
rav-tracing  argument,  Gurnett  et  al.  [1986]  suggested  that  the  broadband  emissions  are 
whistler  mode  radiation  from  the  electron  beam. 

Several  studies  have  attempted  to  explain  how  to  produce  the  broadband  whistle  radi¬ 
ation  by  a  several  keV  electron  beam.  By  calculating  the  radiated  power,  Farrell  et  al. 
[1988]  proposed  that  the  electron  beam  generated  the  radiation  through  a  coherent  process, 
since  the  measured  wave  power  is  several  orders  higher  than  the  power  expected  from  in¬ 
coherent  Cerenkov  radiation.  They  further  suggested  that  electron  bunches  created  by  an 
electrostatic  beam-plasma  instability  radiate  the  whistle  mode  noise.  However,  they  have 
not  shown  how  the  electron  beam  can  excite  the  electrostatic  beam-plasma  instability  below 
the  cyclotron  frequency. 

Solving  the  electrostatic  dispersion  equation  of  a  finite-size  beam  for  parameters  relevant 
to  the  Spacelab  experiments,  we  previously  found  that  the  beam-plasma  instability  can 
generate  broadband  electrostatic  beam  mode  with  frequencies  below  the  election  cyclotron 
frequency  [Lin  and  Wong,  1986).  The  maximum  growth  rate  of  the  beam  mode  in  this 
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frequency  range  is  about  10-'1u>pf.  To  explain  whistler  wave  emissions,  we  suggested  that 
the  broadband  electrostatic  waves  convert  to  whistler  waves  at  the  electron  beam  edge.  Re¬ 
cently,  a  computer  simulation  using  an  electromagnetic  particle  code  showed  that  a  spatially 
confined  thin  electron  beam  can  directly  excite  whistle  mode  waves  near  the  resonance  cone 
through  Landau  resonance  [Omura  and  Matsumoto,  1988].  This  simulation  also  suggested 
that  bunched  electrons  could  not.  radiate  oblique  whistler  waves  because  large  amplitude 
electrostatic  waves  destroy  the  bunching  in  a  short  time  during  the  quasilinear  phase. 

The  objective  of  this  study  is  to  investigate  the  various  wave  modes  that  can  be  excited 
by  a  finite-radius  electron  beam,  with  an  emphasis  on  wave  frequencies  below  the  electron 
cyclotron  frequency.  To  achieve  this  objective,  we  solved  the  electromagnetic  dispersion 
equation  of  a  finite-radius  beam  immersed  in  a  uniform  background  plasma.  For  compari¬ 
son  purposes,  we  also  solved  the  dispersion  equation  of  a  homogeneous  beam.  The  numerical 
solutions  indicate  that  a  keV  electron  beam  can  drive  the  beam  and  whistler  modes  unsta¬ 
ble,  regardless  of  whether  the  beam  is  homogeneous  or  has  a  finite- radius.  Although  the 
beam  mode  is  excited  for  frequencies  up  to  the  electron  pHsma  frequency,  only  waves  with 
frequencies  below  the  electron  cyclotron  frequency  can  propagate  outside  the  beam.  The 
parallel  wavenumbers  of  the  whistler  waves  excited  by  the  finite-radius  and  the  homogenous 
electron  beams  are  similar,  suggesting  that  the  finite-radius  electron  beam  also  excites  the 
whistler  waves  near  the  resonance  cone.  These  results  explain  that  a  very  narrow  electron 
beam,  such  as  those  of  the  Spacelab  experiments,  can  still  excite  whistler  waves  in  the 
resonance  cone. 
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III.  Formulation 


A.  Basic  Equations 


This  section  describes  the  derivation  of  the  dispersion  equation  of  a  finite-radius  electron 
beam  in  a  uniform  cold  background  plasma.  We  assume  that  the  ambient  magnetic  field 
B0  =  B0z  is  in  the  ^-direction  and  that  the  electron  beam  propagates  along  the  magnetic 
field. 

We  consider  an  electron  beam  with  a  density  74  and  a  radius  r„  centered  at  the  origin 
r  =  0.  We  assume  the  background  electrons  have  a  density  n0  outside  the  beam  and  nc  — 
(n0  —  Tib)  inside  the  beam  so  that  the  background  ions  distribute  uniformly,  as  illustrated 
in  Figure  1.  The  relevant  equations  are  the  fluid  equations 


drij 

~dt 

£3. 

dt 

and  the  Maxwell  equations 


4-  V  ■  (rijVj)  =  0 

+  Vj  VV}  =  --(£  +  V,  X 
m 


IdB 

V  x  £  = 

c  dt 

1  dE  477 

V  x  B  = 

c  dt  c  j 

<1 

II 

477  £  tjTlj 
j 

V  •  B  = 

0 

(1) 

B0  +  B))  (2) 

(3) 

i  Vj  (4) 

(5) 

(6) 


where  rij  and  Vj  are  the  number  density  and  velocity  of  species  j,  respectively.  Here  the 
index  j  =  b,  p  refers  to  the  beam  and  background  electrons,  respectively.  For  the  background 
electrons,  vop  =  0. 

The  electromagnetic  dispersion  equation  of  a  finite-radius  beam  has  been  derived  by 
Shoucri  and  Kitsenko  [1968].  Here,  we  briefly  outline  the  procedure  for  obtaining  the 
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dispersion  equation.  To  linearize  (l)-(6),  we  assume 

=  vojz  -I-  Si T} 

B  =  B0z  +  bB 
E  =  6E 
nj  =  rtj(r) 

where  the  zeroth  order  velocity  Vj  is  constant  and  in  the  magnetic  field  direction.  For  the 
background  electrons,  Vop  =  0.  In  the  derivation  given  below,  we  define  f le  =  (|e| B0)/{mec), 
and  u'pj  =  (47rrjje2)/me.  Using  the  cylindrical  coordinates,  we  represent  the  perturbed 
quantities  bf  as 


bf  =  6f(r)e‘<'k^~wt-m0) 


where  k (|  is  the  parallel  wave  number,  8  is  the  azimuthal  angle,  and  m  is  an  integer.  Equa¬ 
tions  (1)— (2)  then  become 


i(uj  —  fc||i>0j  )6rij  = 

n,  d  .  .  ,  -  .71771,  .  . 

O  (rbVjr )  4-  tkunjbVj.  4-  ?  bvjg 

r  or  r 

,  .  drij 

+  Sv‘r~fr 

(T) 

II 

u 

P 

tv-o 

J? 

p 

1 

•  Pi 

1 

e 

m 

bEr  -  - VojbBo  +  -bv'jgBg 
c  c 

(8) 

-i(u>  -  k\\Vo})6vj$  = 

e. 

rn 

6Ee  +  -v0j6Br dujr  B0 

c  c 

(9) 

-i(u  -  k\\v0J)6v'jt  = 

~6EZ 

m 

(10) 

From  (3),  we  obtain 


-6E,  -  ill  BE, 

(11) 

ujr  u 

^6Er  +  -~bEt 
u  u?  err 

(12) 

6Br 

c 

6Bg 


Substituting  (11)  and  (12)  into  (7)- (10),  we  derive  the  perturbed  velocity  quantities  as 

bl’)r  =  RJTbET  +  RjgbEff  4  RJzbEz  (13) 
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6v}$  =  Qjr6Er  +  ®}gt>E$  -f-  QJzbE;  (14) 

Sv~s  =  ZJtSEt  (15) 


where 


Rjr 

it  (u;  -  liiUoj)2 

mu:  (u:  -  !||t>0j  )2  -  ft2 

R}6 

eil  (lu  -  k\\v0]) 
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Rjz 

e  i’oj 
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Qjr 

= 

-Rje 
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Rjr 

eJt 
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if  1 
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— 
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d  mil 

dr  r 


d  m 

-Q— +  (u;- 


The  perturbed  current  density  6  J  is  defined  as 

=  ^2e(6njVOJz  +  TijSvj) 


=  a  -6E 


(16) 

(17) 


where  a  is  the  conductivity  tensor.  The  dielectric  tensor  is 


r  4^r 

*tj  —  oXJ 

tuJ 


(18) 


Substituting  (7),  ( 13)— ( 15)  into  (18).  we  get  the  expression  for  the  dielectric  components 

err  =  1  _  £  ^£i _ (^  ~  k\\v°^2 _ 


V  (u 

—  _  V'  lU>Pl  il{u>  —  l'||t'oj  )2 
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f~n  = 
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dr 
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fee 
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-Er0 
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mil  1  d 


-uJ 
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B.  Normal  Components 


Suhl  and  Walker  ( 1954)  had  expressed  the  system  of  equations  in  terms  of  the  tangential 
compounds  of  the  electric  and  magnetic  fields  8EZ  and  8BZ.  For  the  axially  symmetric 
plasma  oscillations  in  a  cylindrical  geometry,  we  express  below  the  normal  components  of 
8E  and  8B  through  bE:  and  8BZ. 

We  first  linearize  (4)  and  (6)  to  obtain 


l^6Bt  -  ik\\8Bg  =  — t  —  ^2  erj8Ej 

) 

(19) 

>*||W  -  =  ->^^e,,6E, 

(20) 

-%-(r6B,)  -  1 ™Br  =  -,--£'',6E, 
r  or  r 

j 

(21) 
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(22) 
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Substituting  (11)  and  (12)  into  (19)  and  (20),  we  get 


l  -l~c*r8  ^  \  (  6Er  \ 

i  -*7 t6r  +  *“ir  j  V  6Eo  ) 

=  (  <?*.  +  ££  T  ^ 

V*>.  +  ^  -4/ 

Solving  (23),  we  derive 
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Given  6Er  and  8E$,  6Br  and  8Bg  are  deduced  from  (11)  and  (12) 


(23) 


To  make  the  derivation  of  the  system  of  equations  easier,  we  rewrite  the  expressions  of 
8Er  and  8E$  using  the  following  definitions: 

_  =  1  v  -  A'n^)2 
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From  these  definitions  we  get 
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where  JVy  =  fcy  cj^j.  The  expressions  for  6Er  and  6E$  then  become 
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We  can  rewrite  SB±  using  (26),  as 


6Bx  —  -^||(el|  *  x 
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Equations  (25)  and  (26)  can  be  combined  into 
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IV.  Dispersion  Equation 


A.  System  of  Equations 

We  now  proceed  to  derive  the  equations  from  the  tangential  components  of  Maxwell 
Equations 


SB.  =  — — Vx  x  6E± 

u > 


1C 


— Vx  x  SB±  =  SDZ 


=  eZTSET  +  tzB&Ee  ezz6Ez 


(29) 


(30) 


Taking  the  curl  of  (29),  we  deduce 


i.  ? 
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(31) 


,  alt.-  A  1  d  d  m2 

where  Ax  is  the  Laplacian  operator  Ax  =  --r-r— - r  • 

r  Or  Or  rl 

Taking  the  curl  of  (30)  and  using  (28),  we  deduce  (30)  to 


[— (A'|f  4-  &)(A  +  1)  +  C2]Ax6-E?z  +  i(  A|^  4-  b)£A±6Bz 


-  I-3 

—  All 
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(32) 


After  rewriting  (31)  and  (32),  the  system  of  equations  is 


A  ±6EZ  +  K6EZ  +  C6BZ  =  0 
A  x6Bz  +  V6Ez  +  Q6Bz  =  0 


(33) 

(34) 


where 
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V 


Q 


.  k?.  Y 

*TT -z(a  +  b+  Cl  ) 

N\\z 

j2(a  -  Nj)Z  +fY  +  (b  +  Nj)d(f  -d)  +  c2cP 


and 

Y  =  (6  +  N\\)d  +  (a  —  N^)f 
Z  =  (fc  +  iVy  )(a  fe)  —  c2(c  —  JVj|  ) 

We  next  introduce  two  new  variables  'I'i  and  $2 

=  6EZ  +  ihxA6Bz 
^2  =  8EZ  -j-  i\.2A8B  z 

where  Ai  and  A2  are  the  roots  of  the  quadratic  equation 

A2  +  FA  -  1  =  0 

The  coefficient  F  is  given  as 

F  =  {JVjf(c,  -c2)-ATf[a(Cl  +  fe-2C2)  +  62  +  ((/-/)2] 
-  b(d  -  f  )2  +  (a  +  b)(ab  +  /2)  -  c2(a2  -  d2)}/  N\\ 
/y[(a  +  b  +  c)[(b  +  N*)d  +  f(a  -  A/p] 

The  coefficient  >4  is  determined  from  the  condition  that 

A2V  =  -C 


Hence 


A  = 


1 

\/a  +  6  +  ci 


(35) 

(36) 


(37) 


(38) 


(39) 
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Multiplying  (33)  by  i\A  and  adding  it  to  (34),  we  obtain 

Vi  +Tf*i  =  0 

Vx^2  +  T22^2  =  0 

with 

T2  =  fC  +  i^AV 
Tj  =  1C  -f-  iAjA'P. 


(40) 

(41) 


(42) 

(43) 


Using  (37),  (42)  and  (43),  one  can  show  that  T  satisifies  the  equation 

_  T’${-N!(2a  +  26  +  c,  +  *) 
+  iV2[a(2a  +  b  +  Ci  +  2 c2)  -  b2  -  (d  -  /)4]  +  (a  +  b)(ab  +  /2) 

-  b(d  -  /)2  -  c2(a2  -  d2)}  +  tf04(a  +  6+  c,)[(a  -  N2)  -  d2]  =  0 


(44) 


where  K0  =  uj/c.  In  the  case  of  an  unbounded  plasma  and  beam,  T  corresponds  to  a  normal 
component  of  the  wave  vector  Kx  and  to  the  dispersion  equation  for  an  homogenous  plasma 
and  beam. 

The  components  6EZ  and  6BX  are  related  to  the  functions  'I'1  and  2  by 


6Et 

8BZ 


A 2$h  -  Ai^2 
A2  —  Ai 
*  -  4>2 

/lAj-f, 


(45) 

(46) 


For  axially  symmetric  oscillations,  a  solution  of  (40)  and  (41)  in  general  has  the  form 


^*1.2  =  ■Sj>2Jro(7i,2r)  -f  Gi,2?o(Ti>2r)  (47) 

where  J0  and  F0  are  the  Bessel  and  Neumann  functions,  respectively.  Si,2  and  Gji2  are 
constants. 
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B.  Boundary  Conditions 


The  boundary  conditions  are  obtained  by  integrating  Maxwell  equations  across  the  beam- 
plasma  interface  from  ra  —  e  to  ra  +  t  with  f  — >0.  The  boundary  conditions  are 

(6EZ)  =  0 
(6  Eg)  =  0 
(SB,)  =  0 

(SB,)  =  (-jj-SEr)  -  (—j^-SB,)  (48) 

_ _ 

o/(o;  -  fc||Uo;)[(w  -  h Voj)2  - 

_ l_/ST' _ Ulj^k\\Vlj _ 

N\\  \  i  W2(W  ~  V<v)l(W  -  k\\Voj)2  “ 

The  brackets  {)  denote  the  difference  in  the  values  across  the  beam  surface.  For  example, 

(SE,)  =  SEf  |r=r.+0  -  6El\r=ra- o 


The  superscripts  O  and  I  mean  the  quantity  is  evaluated  outside  and  inside  the  beam, 
respectively.  The  boundary  conditions  indicate  that  SEt,  SEg  and  SB,  are  continuous 
at  the  beam  surface  but  that  SBg  is  not.  Using  (45)  and  (46),  we  write  the  boundary 
conditions  in  terms  of  the  functions  and  'I' 2-  The  first  two  boundary  conditions  mean 
that  the  tangential  components  of  the  electric  fields  are  continuous  at  the  interface. 


A2*i  -  Ai»a. 

'  Aj  -  A,  ' 
,  k0  _  i  d*A 

'(A2-Ai) A\T?  dr  Tl  dr  )] 


MAj-A,' 
.  iK0  /A ,5*,  _  Aa  0*2 V 
\A2-Ai)A*\T?  dr  T2  dr  ) 


0 

0 

0 

0 


(49) 
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Parameters  are  different  inside  and  outside  the  beam  because  density  and  beam  velocity 
are  different  (see  Figure  1). 

C.  Dispersion  Equation 

Because  electric  and  magnetic  fields  should  be  finite  on  the  Z  axis,  the  functions 
inside  the  beam  should  be 


<2  =  S'u2Jm(Tt2r).  (50) 

The  fields  outside  the  beam  should  be  free  propagating  waves  at  a  large  distance  (r  — >  oo): 
therefore,  the  functions  2  outside  the  beam  are 


* 


o 

1,2 


(51) 


where  H^(T°2r)  is  the  mth  order  Hankel  function  of  the  first  kind.  S[2  and  G°2  are 
constants. 

The  boundary  conditions  give  four  equations  for  the  four  constants  S[2  and  G®2.  Requir¬ 
ing  the  determinant  of  these  equations  to  be  zero,  we  obtain  the  dispersion  equation 


where 


Dn 

Du 

Du 

Du 

D2\ 

D22 

D23 

D24 

Dz\ 

d$2 

D33 

D34 

Z?4i 

Du 

D*  3 

Du 

(52) 


Du  =  A  [MTlr) 
Du  =  A  {JmiT'r) 
Du  =  A  MWr) 
Du  =  A?fJ<|>(7?r) 
Dn  =  Jm{T[r)IA 1 
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D„  =  Jm{T{r)/A‘ 

b*>  =  HllHifr  )/A° 

Oj.  =  ffSWr)/*0 

I>3,  =  qtfrjTlr)  -  -p'7m(T,V) 

r 

On  =  llT;j'm(T]r)  -  jp{jm(T{r) 

O33  = 

r 

D„  =  pjT/j'JT/r)  +  -w'jm(r'r) 

r 

0«  =  nWJ'JTir)  +  ™w[j„(T‘r) 
0,3  =  90tfHW'(T?  r)  + 

D„  =  +  —  u.?/ft;>(7?r) 

r 


The  dispersion  equation  (52)  contains  the  variables  p,,  90  90  and  which  are  defined  as 

[(a  -  W,f)(6  +  JVJ)  +  /dJArr  -  JV„d/>l 
Pl’2  ~  (a  —  iV,?)2  —  d2  r_ 


91,2 


cAT^ 


VAu 
51,2  ”  cA2^2, 

(6  +  iV,f)2  +  2/(6  +  JV,?)(a  -  JVff )  +  pd  -  l\(a  -  N||2)2  -  d2]  Au 
“'W  "  (a  -  iV,2)2  -  rf2  *n 

l(q-^)(t  +  iV{)  +  /d 
"A  (a-iVjf)2-^ 


where 

,  _ 

J  W2(u;  -  fc|jv0j)[(u;  -  fc||Voj)2  -  ft2] 

In  general,  the  variables  depend  on  n}  and  v0}  and  thus  vary  from  inside  the  beam  to 
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outside  the  beam.  We  thus  use  the  superscripts  I  and  0  to  distinguish  their  values  inside 
and  outside  the  beam,  respectively,  in  Section  II. 
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V.  Description  of  Numerical  Code 

This  section  describes  the  numerical  code  EMFBP  that  solves  the  dispersion  equation 
(52).  EMFBP  contains  the  main  program  and  the  five  subroutines  FCN,  FSBELEC,  JHC, 
ROOT  and  DFCN.  Although  the  dispersion  equation  is  quite  complicated  and  requires 
double  precision  to  evaluate,  the  structure  of  this  program  is  simple.  The  program  computes 
the  complex  roots  of  the  dispersion  equation,  which  is  a  complex  algebra  equation,  for  given 
input  parameters.  The  root-solving  routine  iterates  until  the  error  in  solving  the  dispersion 
equation  is  less  than  a  prescribed  value,  10-4.  The  program  then  outputs  the  root  for  each 
k  value. 

A.  Main  Program 

The  main  program  primarily  reads  the  input  variables,  which  are  listed  below: 


Table  1.  List  of  Input  Variables 

variable 

name 

definition 

typical 

value 

WR2 

u>pe/fi 

3 

DR 

nb/nc 

0.1 

VB1 

V?P/C 

0 

VB2 

vipt'c 

0 

VBB 

VL/C 

0 

RBA 

Mpe^a/ Vob 

0.111 

KI 

fc||C/fl 

40 

DK 

Afc||C/fl 

1 

NK 

number  of  k  mode 

10 

GR 

real  (Z) 

0.33 

GI 

Imaginary  (Z) 

0.05 

The  variable  WR2  is  the  ratio  of  the  electron  plasma  frequency  to  the  electron  gvrofrequency. 
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As  shown  in  the  density  model,  the  variable  DR  is  the  ratio  of  the  beam  density  to  the 
background  density  outside  the  beam.  The  variables  VBl  and  VB2  are  the  drift  velocities 
of  the  background  electrons  outside  and  inside  the  beam,  respectively.  In  this  study,  we 
have  assumed  VBl  and  VB2  to  be  zero.  The  variable  VBB  is  the  beam  velocity  in  the  unit 
of  speed  of  light.  For  a  1  keV  electron  beam,  VBB  is  0.01.  The  beam  radius  defines  the 
normalized  variable  RBA;  for  a  beam  radius  of  2  m  and  a  beam  energy  of  1  keV,  RBA  is 
0.0111.  The  main  program  also  reads  the  initial  value  of  fc||c/Q  for  solving  the  dispersion 
equation.  DK  is  the  increment  in  fcjjc/Q  and  NK  is  the  number  of  times  the  program  solves 
the  dispersion  equation.  The  program  requires  an  initial  guess  of  the  root  Z,  where  Z  is 
ui/Q.  The  input  variables  are  contained  in  a  common  block  /MON/  for  passing  arguments 
to  other  subroutines. 

After  reading  the  input  variables,  the  main  program  calls  the  subroutine  ROOTF,  which 
uses  Muller’s  method  for  solving  a  complex  algebra  equation.  In  ROOTF,  the  subroutine 
calls  DFCN  for  finding  the  roots  and  FCN  for  computing  the  dispersion  equation  (52). 

B.  Computation  of  the  Dispersion  Equation 

Complex  function  FCN*16(Z)  function  computes  the  dispersion  equation  (52)  for  a  given 
Z  and  returns  the  result  in  a  complex  value.  FCN  calls  two  subroutines  FSBELEC  and 
JHC.  The  subroutine  FSBELEC  computes  the  electron  response  terms  a,  b ,  cj,  c2,  d,  and 
/,  given  in  Section  III  A,  and  l ,  given  in  Section  IV  C. 

The  subroutine  JHC  computes  the  mth  order  Bessel  function  Jm(Z)  and  Hankel  function 
of  the  first  kind  H^}(Z).  The  subroutine  JHC  has  input  arguments  Z  and  m,  where  Z  is 
complex  and  m  is  integer.  The  subroutine  computes  Jm(Z )  and  H^(Z)  according  to  the 
expressions  given  in  Handbook  of  Mathematical  Functions  [Abramowitz  and  Segun,  1968]. 
For  the  absolute  value  of  Z  less  than  15,  the  routine  uses  the  series  expansion  to  compute 
Jm(Z)  and  H£\Z).  For  the  absolute  value  of  Z  greater  than  15,  the  routine  uses  Hankel’s 
asymptotic  expansion. 


VI.  Numerical  Results 


This  section  describes  the  numerical  solutions  of  the  dispersion  equation  for  a  homoge¬ 
neous  beam  and  for  a  finite-radius  beam.  We  assumed  that  the  cold  electron  beam  has 
a  beam  energy  of  1  keV  and  a  density  1%  of  the  ambient  density,  nt,/(n,.  +  a*,)  =  0.01. 
This  assumption  was  based  on  the  Spacelab  2  experiment  parameters,  which  had  a  1-keV. 
50-mA  electron  beam  [Gurnett  et  al.,  1986].  The  ambient  plasma  density  nc  at  the  Shuttle 
altitude  is  generally  near  105  cm-3,  corresponding  to  a  plasma  frequency  of  2.85  MHz.  Since 
the  Spacelab  2  measured  the  electron  plasma  waves  near  3  MHz,  we  assumed  the  electron 
plasma  frequency  to  be  3  MHz  in  solving  the  dispersion  equations.  We  also  used  we/17e  =  3 
because  the  electron  cyclotron  frequency  is  about  1  MHz.  In  solving  the  dispersion  equation 
of  a  finite-radius  beam,  we  assumed  the  beam  radius  to  be  2.5  m,  according  to  Gurnett  et 
al.  [1986]. 

Figs.  2  and  3  give  the  solutions  of  the  dispersion  equation  for  the  homogeneous  beam. 
Fig.  2  shows  the  wave  frequency  ujr  as  a  function  of  kc/Qt  for  various  wave  normal  angles. 
The  beam  plasma  mode  has  frequencies  linearly  proportional  to  kcj 17e.  intersecting  the 
plasma  frequency  at  uje/Sle  =  3,  as  indicated  on  the  vertical  axis.  The  dashed  lines  in 
Fig.  2  represent  the  whistle  mode,  which  has  frequencies  approaching  Qe  cos  8.  Fig.  3  gives 
the  growth  rate  u;,  for  both  the  beam  plasma  and  whistler  modes.  As  expected,  the  beam 
plasma  mode  has  a  large  growth  rate,  which  decreases  with  angle  (solid  lines).  For  example, 
the  maximum  value  of  u\/17e  is  0.4  for  0  =  0°  and  0.25  for  6  =  60°.  Note  that  the  growth 
rate  of  the  beam  plasma  mode  at  large  angles  has  dropped  sharply  at  certain  kc/?lt.  For 
example,  u\/S7e  for  0  =  60°  drops  below  10-2  for  kc/Ue  between  20  and  32,  corresponding 
to  wave  frequencies  near  the  electron  cyclotron  frequency.  The  decrease  in  growth  rate  is 
due  to  the  cyclotron  resonance  interactions  of  the  electron  beam.  Fig.  3  also  shows  that  the 
cold  electron  beam  can  excite  whistler  waves  at  large  normal  angles  with  smaller  growth 
rates  (dashed  lines).  FYom  Figs.  2  and  3,  we  note  that  the  frequency  of  unstable  whistler 
waves  occurs  at  ft*  cos  9  where  the  frequency  is  insensitive  to  k,  suggesting  that  the  electron 
beam  excites  the  whistler  mode  at  the  resonance  cone.  Fig.  3  indicates  that  the  growth  rate 
of  the  whistler  mode  is  about  0.025ft,  for  B  =  60°  and  0.01217,  for  6  =  30°.  For  0  =  60°,  this 
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corresponds  to  an  e-folding  time  of  approximately  6  p s.  Although  the  temporal  growth  rate 
of  the  whistler  wave  is  small  in  comparison  with  the  beam  mode,  the  spatial  amplification 
of  the  whistler  wave  can  be  significant  because  of  its  small  group  velocity  at  the  resonance 
cone. 

Figs.  4  and  5  give  the  solutions  of  the  dispersion  equation  of  the  finite-radius  beam.  In 
solving  (52),  we  restrict  our  attention  to  the  axially  symmetric  case  (i.e.,  m  =  0)  only. 
Since  we  are  mainly  interested  in  waves  with  frequencies  below  fle,  we  have  neglected  the 
cvclotronic  mode,  which  has  been  discussed  previously  by  Le  Queau  et  al.  [1981].  Because 
the  perpendicular  wavenumber  is  no  longer  a  constant,  these  figures  plot  the  frequency 
and  growth  rate  versus  the  normalized  parallel  wavenumber  k\\c/Q,e.  Like  the  homogeneous 
beam  case,  the  finite-radius  electron  beam  is  unstable  to  the  beam  plasma  mode  with 
wave  frequency  ur  =  A.||t’0«>  (Fig.  4).  The  growth  rate  is  also  large  when  k^c/£le  >  20, 
which  corresponds  to  u-v/fl,  >  1.  However,  except  near  the  electron  plasma  frequency,  the 
solutions  indicate  that  the  argument  of  the  Hankel  function  for  the  beam  mode  has  its 
imaginary  part  larger  than  the  real  part  (j/m(7^r)|  >  )( )  when  u}r/Qe  >  1  (Fig. 

4).  This  means  that  the  waves  with  frequencies  <  uj  <  u>r  are  radially  confined  inside 
the  beam.  The  reason  for  the  confinement  is  that  the  ambient  plasma  outside  the  beam 
does  not  support  a  normal  mode  for  Qe  <  u  <  we.  For  u>  <  Qe,  the  beam  mode  has  a 
smaller  growth  rate  about  0  01flr,  but  the  argument  of  the  Hankel  function  has  a  larger 
real  part  than  imaginary  part.  Therefore,  the  beam  mode  excited  with  frequencies  below 
Qe  can  propagate  away  from  the  beam  because  the  ambient  plasma  has  the  whistler  mode 
as  its  normal  mode.  The  drop  in  the  growth  rate  at  the  cyclotron  frequency  is  due  to  the 
cyclotron  resonance  interactions  of  the  electron  beam,  similar  to  the  homogeneous  beam 
case. 

Besides  the  beam  mode,  (52)  also  has  solutions  for  the  whistler  mode  shown  in  Fig.  5. 
This  figure  indicates  that  the  frequency  approaches  for  large  k^c/Qr-  The  solutions  for 
ojr/Q.(  <  0.02  are  near  the  lower  hybrid  frequency  and  thus  are  not  accurate  because  we  have 
neglected  the  ion  contributions  in  (52).  Although  the  real  frequency  appears  as  a  single 
curve,  we  have  actually  obtained  several  unstable  branches,  all  with  very  close  frequencies 
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but  different  growth  rates.  Fig.  2  indicates  that  the  whistler  mode  for  the  homogeneous 
beam  has  a  frequency  for  each  angle.  Therefore,  at  a  given  parallel  wavenumber,  the  whistler 
mode  has  a  continuous  frequency  spectrum  as  the  perpendicular  wavenumber  varies.  For 
a  finite-radius  beam,  the  perpendicular  wavelength  inside  the  beam  is  quantized  and  thus 
the  whistler  mode  has  a  set  of  discrete  frequencies  as  shown  in  Fig.  5.  We  have  plotted 
two  branches  with  the  largest  growth  rates,  which  are  smaller  than  the  growth  rates  of  the 
beam  mode. 

For  the  finite-radius  beam  case,  one  cannot  easily  determine  whether  the  electron  beam 
excites  whistler  waves  at  the  resonance  cone,  since  the  perpendicular  normal  vector  is  not 
constant.  However,  we  note  that  fc||c/f!e  is  the  same  for  the  finite-radius  and  homogeneous 
beams  cases.  Fig.  2  indicates  that  the  homogeneous  beam  excites  whistler  waves  at  k\\c/Qe  = 
24  for  kc/£le  =  48  and  6  —  60°,  whereas  the  finite-radius  beam  excites  whistler  waves  for 
*||c/fle  ranging  from  18  to  25,  according  to  Fig.  5.  The  parallel  wavelength,  therefore, 
suggests  that  the  finite-radius  beam  also  excites  whistler  waves  near  the  resonance  cone. 
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VII.  Summary 


Our  numerical  calculations  indicate  that  the  keV  electron  beam  injected  from  the  Shuttle 
with  a  narrow  radius  can  directly  excite  whistler  waves  propagating  at  the  resonance  cone. 
This  result  thus  supports  the  work  of  Omura  and  Matsumoto  [1988],  which  showed  the 
direct  excitation  of  whistler  waves  by  using  computer  simulations.  Although  the  beam 
mode  is  excited  for  frequencies  up  to  the  electron  plasma  frequency,  only  the  waves  below 
the  electron  cyclotron  frequency  can  propagate  outside  the  beam.  As  suggested  by  Farrell  et 
al.  [19S8],  the  electrostatic  waves  of  the  beam  mode  at  the  whistler  frequencies  could  bunch 
electrons,  which  then  radiate  whistler  waves.  However,  one  cannot  rule  out  the  possibility 
of  mode  conversion  at  the  beam  edge,  as  proposed  by  Lin  and  Wong  [1986].  Electrostatic 
waves  of  the  beam  mode  could  convert  into  whistler  waves  at  the  beam  edge  due  to  density 
discontinuity  and  propagate  away  from  the  beam.  Future  studies  are  required  to  resolve 
these  questions. 
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Figure  Captions 


Figure  1.  Physical  configuration  of  a  finite-radius  electron  beam  immersed  in  a  uniform 
plasma. 

Figure  2.  Frequency  of  the  beam  mode  (solid  curve)  and  the  whistler  mode  (dashed  curve) 
at  a  given  wave  normal  angle  8  as  a  function  of  kc/Qe. 

Figure  3.  Growth  rate  of  the  beam  mode  (solid  curve)  and  the  whistler  mode  (dashed 
curve)  at  a  given  wave  normal  angle  6. 

Figure  4.  Frequency  and  growth  rate  of  the  beam  mode  for  a  finite-radius  electron  beam  as 
afunction  of  k\\c/Qe.  The  parameters  are  nt,/nc  =  0.1,  Eb  =  1  keV,  ra  =  2  m,  andu?e/Qe  =  3. 

Figure  5.  Frequency  and  growth  rate  of  the  whistle  mode  for  a  finite-radius  electron  beam. 
The  parameters  are  the  same  as  those  for  Figure  4. 
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Abstract 

The  full  electromagnetic  dispersion  equation  of  a  homogeneous  beam  and  the 
electrostatic  dispersion  equation  of  a  finite-size  electron  beam  are  solved  for 
parameters  relevant  to  DC  electron  beam  experiments  conducted  from  the 
Space  Shuttle.  The  numerical  solutions  of  the  homogeneous  dispersion  equa¬ 
tion  indicate  that  the  electron  beam  mode  with  frequency  u>  =  has  large 
growth  rates  over  a  wide  range  of  frequencies  and  angles,  whereas  whistler 
waves  have  small  growth  rates  at  frequencies  near  the  electron  cyclotron 
frequency  fle.  At  oblique  propagation  and  at  low  frequencies  (w  <  0.1fte), 
the  beam  mode  is  coupled  to  whistler  waves  and  has  a  significant  electro¬ 
magnetic  component.  Because  of  the  finite  6ize  effects,  the  electron  beam 
instability  has  a  minimum  growth  rate  at  the  electron  cyclotron  frequency. 
Our  results  do  not  favor  direct  excitation  of  whistler  waves  in  the  ionosphere 
through  linear  instabilities  of  keV  electron  beamfe.  To  explain  whistler  wave 
emissions  during  the  Shuttle  electron  beam  experiments,  it  is  suggested  that 
the  electron  beam  first  excites  broadband  electrostatic  waves,  which  are  sub¬ 
sequently  converted  to  whistler  waves  at  the  beam  edge. 
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1  Introduction 


Different  types  of  waves  have  been  observed  in  association  with  the  DC  elec¬ 
tron  beam  emission  experiments  in  the  Spacelab  1  and  Spacelab  2  missions. 
During  the  Spacelab  1  mission,  broadband  emissions  at  frequencies  less  than 
the  electron  cyclotron  frequency  {u>  <  Qe)  with  both  electric  and  magnetic 
components  were  detected  by  wave  instruments  on  board  the  Shuttle  during 
the  firing  of  a  5  keV  and  300  mA  electron  beam  (Taylor  et  al.,  1985).  Dur¬ 
ing  the  recent  Spacelab  2  mission  the  Plasma  Diagnostics  Package  (PDP) 
performed  a  fly-around  of  the  Shuttle  at  distances  of  up  to  300  meters  while 
a  DC  electron  beam  with  a  beam  energy  of  1  keV  and  current  of  50  mA 
was  being  emitted  from  the  Shuttle.  The  PDP  detected  broadband  emis¬ 
sions  at  u>  <  fle  with  both  electric  and  magnetic  components  and  a  narrow 
band  electrostatic  emission  near  the  plasma  frequency  wpe  (Gurnett  et  al., 
1986).  The  broadband  emissions  had  a  clear  funnel  shape  in  the  frequency¬ 
time  spectrogram,  similar  to  the  observations  of  auroral  hiss  (Gurnett  et  al., 
1983).  The  broadband  emission  is  believed  to  be  caused  by  whistler-mode 
emission  from  the  beam  propagating  near  the  resonance  cone. 

Electron  beams  are  known  to  radiate  whistler  waves  through  sponta¬ 
neous  emissions,  but  the  predicted  intensity  is  generally  several  orders  of 
magnitude  smaller  than  those  observed  (Taylor  and  Shawhan,  1974;  Far¬ 
rell  et  al.,  1986).  We  therefore  examined  the  stimulated  emission  of  waves 
through  instabilities  in  the  whistler  frequency  range.  Using  the  experimental 
parameters,  we  solved  both  the  general  electromagnetic  dispersion  equation 
for  a  homogeneous  beam  and  an  electrostatic  dispersion  equation  for  a  Unite- 
size  beam.  As  expected,  the  solutions  indicate  that  the  beam  mode  is  much 
more  unstable  than  the  whistler  mode.  A  more  surprising  finding  is  that 
the  oblique  beam  mode  has  a  significant  electromagnetic  component  at  low 
frequencies  (w  <  fie).  Based  on  the  instability  analysis,  mechanisms  that 
cause  artificial  electron  beams  to  emit  whistler  waves  in  the  ionosphere  are 
suggested. 
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2  Instabilities  of  a  Homogeneous  Beam 

The  electromagnetic  dispersion  equation  for  the  homogeneous  beam  plasma 
system  is  derived  in  a  representation  based  on  the  electrostatic  and  elec¬ 
tromagnetic  potentials  rather  than  the  electric  field.  The  potential  rep¬ 
resentation  is  useful  for  our  purposes  because  it  simplifies  the  separation 
of  the  electrostatic  and  electromagnetic  perturbations.  Kinetic  effects  are 
included  since  electron  beams  might  be  thermalized  by  beam  plasma  insta¬ 
bilities  after  propagation.  The  dielectric  tensor  used  in  the  calculation  can 
be  transformed  into  the  well  known  formula  of  Stix  (1962)  and  is  therefore 
not  given  here. 

The  ambient  plasma  density  ne  at  the  Shuttle  altitude  is  generally  near 
10scm~3  corresponding  to  a  plasma  frequency  of  about  2.85  MHz.  Since 
the  electron  cyclotron  frequency  is  i.bout  1.0  MHz  in  the  ionosphere,  we 
chose  the  ratio  Vpt/Slc  =  3  in  the  calculation  for  simplicity.  The  two  key 
parameters  in  solving  the  dispersion  equation  are  n&/A\  the  ratio  of  the 
beam  density  nt,  to  the  total  plasma  density  N ,  and  Vt,/a,  the  ratio  of  beam 
drift  velocity  Vj,  to  the  beam  thermal  velocity  a.  The  dispersion  equation  is 
solved  for  the  ratio  of  nt/iV  =  0.1  and  Vt,/a  =  10. 

Figure  1  shows  the  normalized  wave  frequency  u>/fle  and  the  normalized 
growth  rate  y/Cl€  of  the  beam  mode  as  a  function  of  kc/Slt  for  wave  nor¬ 
mal  angle  0  =  0°,  30°  and  60°  (solid  lines).  The  parameters  n^/N  =0.1 
and  Vs/a  =  10  are  used.  The  electron  beam  mode  is  unstable  for  a  wide 
frequency  range  from  zero  to  u>  =  3.4fi„  slightly  above  the  upper  hybrid  res¬ 
onance  frequency  (3.20e).  The  maximum  growth  rate  occurs  at  uj  =  2.6fle 
for  0  =  0°.  Both  the  wave  frequency  and  growth  rate  of  the  beam  mode 
decrease  with  0. 

The  beam  mode  not  only  has  a  large  growth  rate  but  also  has  a  group 
velocity  in  the  direction  near  the  magnetic  field.  For  0  =  30°,  the  angle 
between  the  group  velocity  and  the  magnetic  field  is  found  to  be  only  about 
5°.  Therefore,  the  beam  mode  can  grow  to  large  amplitudes  inside  the 
electron  beam. 

Figure  1  also  shows  the  frequency  and  growth  rate  of  whistler  waves  for 
0  ~  30°  (dotted  line).  The  whistler  wave  frequency  is  close  to  the  beam  mode 
frequency  for  Jtc/fle  <  5.  The  growth  rate  of  whistler  instability  is  much 
smaller  than  the  beam  mode  instability  (note  that  the  whistler  growth  rate 
has  been  multiplied  by  a  factor  of  100  in  Figure  1).  The  instability  occurs 
at  two  frequencies,  near  0.1  fi,  and  0.85  He.  The  growth  rate  at  ut  =  0.lQe 
is  very  small  (only  about  0.0002ftc)  and  not  presentable  in  the  figure.  The 
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Cerenkov  resonant  energy  increases  with  wave  number  k  to  a  maximum  near 
kc/Cle  -  5  and  then  decreases  with  k  after  kc/Cle  >  5.  An  electron  beam 
with  a  beam  energy  of  several  keV  can  resonate  with  whistler  waves  at  two 
wave  vectors  (near  1  and  10  kc/ fie)  and  thus  excite  whistler  waves  at  two 
frequencies. 

To  determine  whether  the  excited  wave  is  electrostatic  or  electromag¬ 
netic,  we  plot  the  ratio  of  the  longitudinal  electric  field  component  El  to 
the  total  electric  field  Ej  (Figure  2).  The  ratio  El/Et  is  1  for  purely 
electrostatic  waves  and  0  for  purely  electromagnetic  waves.  Figure  2  shows 
that  the  beam  mode  is  electrostatic  for  kc/Sle  >  3  and  has  a  significant 
transverse  electric  field  component  ( El/Et  =  0.6  for  kc/£le  =  1).  For 
whistler  waves  at  0  =  30°,  El/Et  is  about  0.4  at  kc/ fle=l  and  about  1  at 
fcc/fle  =10.  Therefore,  the  unstable  whistler  waves  become  electrostatic  for 
kc/ile  >  10.  Since  the  beam  mode  is  generally  believed  to  be  electrostatic, 
it  is  unexpected  that  the  oblique  propagating  beam  mode  has  a  substantial 
transverse  electric  field  component  at  low  frequencies.  The  bottom  panel 
of  Figure  2  indicates  that!  the  dispersion  curves  of  the  beam  and  whistler 
modes  cross  each  other  a i  kc/ile  ss  1,  suggesting  that  the  electron  beam 
mode  couples  to  whistler  graves. 

The  frequency  and  growth  rate  of  the  beam  mode  are  examined  for 
various  values  o  f  n^/N  and  v\>/a.  The  growth  rate  of  the  beam  mode 
decreases  as  the  beam  density  decreases.  At  sufficiently  low  density,  the 
beam  mode  ceases  to  exist  and  the  instability  occurs  only  at  frequencies 
close  to  the  plasma  frequency.  The  effect  of  increasing  the  beam  thermal 
velocity  is  similar  to  that  of  decreasing  the  beam  density. 
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3  Instability  of  a  Finite-size  Electron  Beam 


The  electrostatic  dispersion  equation  of  a  cold  electron  beam  with  a  finite 
radius  injected  along  the  magnetic  field  in  a  cold  uniform  plasma  has  been 
derived  (Le  Queau  et  al.,  1981).  The  electrostatic  dispersion  equation  is 

j,(T)  <°y  /?»>([/) 

W)  -  u  4-  (l) 

where  Jn  is  the  Bessel  function  of  order  n  and  Hn  *  is  the  Hankel  function 
of  the  first  kind  of  order  n.  The  Hankel  function  is  defined  in  terms  of  the 
Bessel  function  of  the  first  kind  Jn  and  second  kind  Yn  as  =  Jn  +  tYn. 
The  arguments  T  and  U  are  defined  as  T  =  k'£a  and  U  =  k^a  .  Here  a  is 
the  beam  radius;  k'£  and  k°fl  are  the  effective  perpendicular  wave  number 
inside  and  outside  the  beam  respectively.  The  following  definitions  are  also 
used: 


In  the  above  equations,  6  is  defined  as  rib/N.  The  equation  of  cj|"  =  0  is  the 
electrostatic  dispersion  equation  for  a  homogeneous  beam  plasma  system. 

Using  the  same  parameters  as  the  homogeneous  beam,  we  solved  Equa¬ 
tion  (1)  numerically.  The  beam  radius  a  is  assumed  to  be  2.5  m,  which  is 
reasonable  for  the  Spacelab  2  experiment  (Gurnett  et  al.,  1986).  Equation 
(1)  contains  solutions  for  three  modes:  Cerenkov,  cyclotronic,  and  surface 
modes  (Le  Queau  et  al.,  1981).  Figure  3  shows  that  the  growth  rate  for 
the  cyclotronic  branch  is  quite  large  (7/fle  *  0.25)  and  the  growth  rate 
of  the  Cerenkov  branch  is  about  one  third  of  the  cyclotronic  branch.  The 
surface  branch  has  the  smallest  growth  rate  among  the  three  branches.  As 
a  comparison,  the  frequency  and  growth  rate  of  the  beam  mode  for  the 
homogeneous  beam  plasma  system  at  0  =  0°  are  also  plotted  in  Figure  3 
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(dotted  lines).  The  frequency  of  the  Cerenkov  branch  is  close  to  the  beam 
mode  frequency.  Note  that  the  finite-size  effects  reduce  the  growth  rate  for 
fcj|Vfc/u>pe  <  1. 

For  a  finite-size  beam  plasma  system,  the  transverse  wave  numbers  are 
quantized.  We  found  that  the  transverse  wave  length  outside  the  beam 
(2 is  about  one  hundred  times  the  beam  radius  for  the  cyclotronic 
branch.  Due  to  lack  of  space,  the  results  will  be  discussed  separately. 

The  growth  rate  at  frequencies  below  f!e  is  very  small  but  interesting. 
Figure  4,  which  presents  the  growth  rate  in  the  logarithm  scale  versus  the 
normalized  frequency  w/wpe,  shows  the  instability  extending  to  frequencies 
below  Qe  =  0.33wpe.  The  growth  rate  has  a  minimum  at  u  =  Qe  and 
increases  with  decreasing  frequency  to  a  local  maximum  around  u>  =  0.25wpe. 
Although  the  normalized  growth  rate  7/wpe  is  only  about  10~3  for  u>  <  De, 
the  wave  amplitude  can  still  grow  to  twenty  times  in  one  milisecond. 
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4  Discussion 


The  full  electromagnetic  dispersion  equation  is  solved  for  parameters  rele¬ 
vant  to  electron  beam  experiments  from  the  Space  Shuttle.  The  numerical 
solutions  indicate  instabilities  for  the  electron  beam  and  whistler  modes. 
The  electron  beam  mode  with  frequency  u  =  k||Vfc  has  large  growth  rates 
over  a  wide  range  of  frequencies  and  angles,  whereas  whistler  waves  have 
small  growth  rates  at  frequencies  near  the  electron  cyclotron  frequency.  At 
oblique  propagation  and  at  low  frequencies  (u  <  0.1$Te),  the  beam  mode 
has  a  significant  electromagnetic  component.  The  dispersion  curve  for  the 
beam  and  whistler  modes  crosses  each  other,  suggesting  that  the  two  modes 
can  be  easily  coupled. 

The  electromagnetic  dispersion  relation  for  a  finite-size  electron  beam  is 
very  difficult  to  solve.  We  have  thus  examined  the  finite-size  beam  effects  in 
the  electrostatic  limit.  The  numerical  solutions  of  the  finite-size  beam  dis¬ 
persion  equation  confirm  the  solutions  of  electromagnetic  dispersion  relation 
of  a  homogeneous  beam,  but  also  exhibit  a  new  and  interesting  feature. 

Our  numerical  solutions  of  the  finite-size  beam  dispersion  equation  agree 
with  the  work  of  Le  Queau  et  al.  (1981)  when  we  use  their  parameters. 
However,  Le  Queau  et  al.  (1981)  solved  the  dispersion  equation  only  for  fre¬ 
quencies  above  the  electron  cyclotron  frequency  fle.  We  obtained  numerical 
solutions  below  the  electron  cyclotron  frequency  and  found  that  the  instabil¬ 
ity  in  fact  extends  to  frequencies  below  fle  with  *  minimum  and  negligible 
growth  rate  at  u  =  fie.  This  result  may  explaiiji  the  Spacelab  2  observa¬ 
tions,  which  show  no  wave  emission  near  fle  (see  Figure  1  of  Gurnett  et  al., 
1986).  Since  the  beam  mode  of  a  homogeneous  beam  has  a  large  growth 
rate  near  u>  =  fle  (see  Figure  3),  the  gap  of  wave  emissions  near  fle  cannot 
be  interpreted  without  invoking  finite-size  effects. 

The  present  study  is  motivated  by  new  experimental  results  from  Space- 
lab  1  and  2  (Taylor  et  al.,  1985;  Gurnett  et  al.,  1986).  In  particular,  the 
remote  sub-satellite  observations  during  the  Spacelab  2  experiment  indi¬ 
cate  broadband  emissions  at  frequencies  less  than  the  electron  cyclotron 
frequency  and  narrow  band  emissions  near  the  electron  plasma  frequency. 
The  broadband  emissions  below  the  electron  cyclotron  frequency  are  sug¬ 
gested  to  be  whistler  waves  propagating  near  the  resonance  cone  (Gurnett 
et  al.,  1986).  Because  the  whistler  mode  instability  occurs  at  oblique  prop¬ 
agation  with  a  small  growth  rate,  it  is  unlikely  that  the  artificial  electron 
beams  fired  from  the  Space  Shuttle  excite  whistler  waves  directly.  However, 
whistler  waves  can  be  generated  by  other  nonlinear  processes  such  as  para- 
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metric  decay  of  the  beam  mode  into  the  whistler  mode.  Beating  of  two  high 
frequency  plasma  waves  can  also  produce  electrostatic  waves  at  frequencies 
below  the  electron  cyclotron  frequency  (La  Queau  et  al. ,  1981).  The  present 
work  suggests  another  mechanism  of  generating  whistler  waves  by  the  di¬ 
rect  coupling  of  the  beam  and  whistler  modes.  Mode  conversion  is  known 
to  occur  easily  when  the  medium  has  a  density  gradient  and  the  dispersion 
curves  of  two  modes  are  crossed  (Stix,  1965).  Therefore,  our  interpretation 
of  whistler  wave  emissions  by  artificial  electron  beams  is  that  the  beam  mode 
at  frequencies  below  the  electron  cyclotron  frequency  is  converted  into  the 
whistler  mode  at  the  edge  of  the  electron  beam. 
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Figure  Captions 


Figure  1  Frequency  and  growth  rate  of  the  electron  beam  for  6  =  0°,30° 
and  60°  as  a  function  of  kc/tte  (solid  lines).  The  frequency  and  growth 
rate  of  whistler  waves  are  also  plotted  as  dotted  curves  for  comparison. 

Figure  2  The  ratio  of  the  longitudinal  to  the  transverse  electric  field  versus 
jfcc/fle  (top  panel).  The  dispersion  curves  of  the  beam  and  whistler 
modes  are  plotted  in  the  bottom  panel. 

Figure  3  Frequency  and  growth  rate  of  plasma  instabilities  of  a  finite-size 
elelctron  beam  as  a  function  of  fc|| The  curves  denoted  by 
1  are  the  Cerenkov  mode.  The  surface  mode  is  represented  by  the 
curves  with  the  numeral  2  and  the  cyclotronic  mode  is  represented 
by  the  curves  with  the  numeral  3.  The  frequency  and  growth  rate  of 
the  electron  beam  mode  for  a  homogeneous  beam  are  also  plotted  as 
dotted  lines. 

Figure  4  Normalized  growth  rate  of  a  finite-size  beam  in  the  logarithm 
scale  as  a  function  of  the  normalized  frequency  u>/u>pe.  The  figure 
shows  that  the  growth  rate  has  a  minimum  at  the  electron  cyclotron 
frequency  and  the  instability  extends  to  frequencies  below  the  electron 
cyclotron  frequency. 
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